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ABSTRACT 

Aims. Spectroscopic, polarimetric, and high spectral resolution interferometric data covering the period 1995-2011 are analyzed to 
document the transition into a new phase of circumstellar disk activity in the classical Be-shell star 48 Lib. The objective is to use this 
broad data set to additionally test disk oscillations as the basic underlying dynamical process. 

Methods. The long-term disk evolution is described using the V/R ratio of the violet and red emission components of Ha and Bry, 
radial velocities and profiles of He i and optical metal shell lines, as well as multi-band B VRI polarimetry. Single-epoch broad-band 
and high-resolution interferometric visibilities and phases are discussed with respect to a classical disk model and the given baseline 
orientations. 

Results. Spectroscopic signatures of disk asymmetries in 48 Lib vanished in the late nineties but recovered some time between 
2004 and 2007, as shown by a new large-amplitude and long-duration V/R cycle. Variations in the radial velocity and line profile 
of conventional shell lines correlate with the V/R behavior They are shared by narrow absorption cores superimposed on otherwise 
seemingly photospheric He i lines, which may form in high-density gas at the inner disk close to the photosphere. Large radial velocity 
variations continued also during the V/R-quiet years, suggesting that V/R may not always be a good indicator of global density waves 
in the disk. The comparison of the polarization after the recovery of the V/R activity shows a slight increase, while the polarization 
angle has been constant for more than 20 years, placing tight limits on any 3-D precession or warping of the disk. The broad H- 
band interferometry gives a disk diameter of (1.72±0.2) mas (equivalent to 15 stellar radii), position angle of the disk (50±9)°and a 
relatively low disk flattening of 1.66±0.3. Within the errors the same disk position angle is derived from polarimetric observations and 
from photocenter shifts across Bry. The high-resolution interferometric visibility and phase profiles show a double or even multiple- 
component structure. A preliminary estimate based on the size of the Bry emitting region indicates a large diameter for the disk (tens 
of stellar radii). Overall, no serious contradiction between the observations and the disk-oscillation model could be construed. 

Key words, stellar interferometry, polarimetry. Stars: circumstellar matter, emission line, Be-Stars; individual: 48 Lib 



1. Introduction 

Circumstellar matter around rapidly rotating classical Be stars 
has been studied since the discovery by Secchi ( 1866[) of e mis- 
sion lines in yCas. Only some 60 years later, |Struve| (| 1 93 1 1) laid 



the foundation to the notion of a flattened structure of this matter. 
It took again more than 60 years until the disk nature was unam- 
biguously proved by the combination of polarimetry and inter- 
ferometry (e.g. Quirrenbach et al. 1997| l. For a general overview 
of classical Be stars and their properties see Porter & Rivinius 
( |2003l l. 

During the last decade, the rate of progress was much ac- 
celerated by the availability for detailed modeling of combined 
high-resolution spatial and spectral observations. This concerns 
especially the disk dynamics. Perhaps the most manifest em- 
pirical indicator of dynamical processes in Be star disks is the 
so-called V/R ratio of the strengths of violet and red emission- 
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* Based partly on observations collected at the European Southern 
Observatory, Chile (Prop. Nos. 383.D-0522, 266.D-5655, 081.C-0475). 



line components. The nature of this variability has long remained 
somewhat enigmatic. Today, it is widely accepted that the long- 
lived and frequently quasi-periodic structural variations of emis- 
sion lines find their basic explanation in one-armed density 
waves of the disk ( Okazaki|1997||Carciofi et al.|2 009 ). Although 
the non-periodicity of the variability makes it difficult to distin- 
guish systematic residuals from ephemeral events, the study of 
density waves and other phenomena triggered by the presence 
of a stellar companion has a strong potential for revealing the 
structure, dynamics, and secular evolution of Be disks (jCarciofi 
[20n] ). 

The activity and even the existence of Be disks typically vary 
on timescales of years, including their temporary absence. The 
gas of the disk is both provided and dissipated by the central 
star, with the disk build-up often resulting from discrete stellar 



mass-loss events/episodes (e.g., Jjaade et al 



2011 Carciofietal. 



20111. It is assumed that the rotationally induced quadrupole 



moment of the stellar gravitational potential is the driving force 
of the disk oscillations (Papaloiz ou et al.||1992 i. In view of ex- 
tended periods without observable V/R activity, this invites the 
question whether variable star-to-disk mass transfer or restruc- 
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Table 1. New and archival 48 Lib spectroscopic datasets in the visual and near-IR range 



Data 
set 


Date 


JD 

2400000+ 


Telescope 


Instrument 


No. of 
spectra 


Resolving 
power 


Spectral 
range [A] 


A 


1995-2003 


49 788-52 725 


ESO 50-cm, Ondi^ejov 2-m 


Heros 


18+18 


20000 


3500-5540, 5800-8600 


B 


2001-07-11 


52101 


ESO/VLT-Kueyen 


UVES 


4 


60000 


3100-10400 


C 


2008-04-03 


54559 


ESOA'LT-Kueyen 


UVES 


9 


60000 


4800-6800 


D 


2007-2011 


54172-55 389 


BeSS database 


LHIRES 


19 


17 000 


6520-6610 


E 


2009, 2011 


55 014,55 695 


OPD/LNA 


ECASS 


4+12 


16000 


3950-4900, 6070-6855 


F 


2009-08-06 


55 049.5445 


Gemini South 


Phoenix 


1 


45 000 


21 600-21 700 



Table 2. High-resolution AMBER and broad-band PIONIER observations of 48 Lib and ^Tau. The elementary integration time 
DIT is given in seconds per frame for AMBER but by scanning frequency in fringe per second for PIONIER. 





Date 


JD 

(+2400000) 


Target 


Instrument 


Baseline 
configuration 


DIT 


Spectral 
band 


Spectral 
resolution 


2011-06-04 
2009-05-13 
2008-11-16 


55 717.1998 
54963.2349 
54754.3453 


48 Lib 
48 Lib 
^Tau 


PIONIER 
AMBER 
AMBER 


DO-GO-Gl-Il 
UT1-UT2-UT4 
UT1-UT2-UT4 


-300 Hz 
3s 
Is 


H 
Bry 
Bry 


15 

12000 
12000 



turings of the disk (e.g., to a more ring-like profile, [Rivinius et"ar 
200 l| l play a role in triggering or quenching density waves. 

Spectro-interferometry has a distinct potential to permit new 
insights to be gathered. After the angular resolution of cir- 
cumstellar disks by broad-band interferometry in the 1990s 



appearance of interferometric visi bihties and phases . Similar HR 
observations of ^ Tau (not used by I 



Stefl et al. 



2009 1 are included 



I not be a singular 



(Quirrenbach et al. 1997 1, the increasing spectral resolution of 
interferometric equipment afforded detailed investigations also 
of the disk dynamics. Gies et al. (2007i determined the size, 
inclination, and orientation of four northern Be-star disks by 
means of broad-band K observations with the CHARA array. 
Later, Meilland et al.| (2007 1 detected an asymmetry in the 



disk of K CMa using the spectro-interferometric capabilities of 
AMBER at medium resolution (MR, R = 1 200). MR AMBER 
observations were also used for a multi-te chnique study of an - 
other bright Be star, ^Tau ( |Stefl et al.|2009||Carciofi et al.|2009 1. 

In the second half of the 20th century, the Ha V/R ratio 
of 48 Lib (HD 142983, B3Ve, V ~ 5.0) went through a numbe r 
of irregular long-term cycles (see, e.g., Hanuschik et al.p995 1. 
Strong, narrow shell absorption lines show that the line of sight 



to demonstrate that this complexity may we " 
peculiarity but rather some commonality only revealed at higher 
spectral resolution. The two stars also have significant long-term 
V/R variations in common. 

Basic geometrical parameters of the disk are derived by 
fitting an axisymmetric Gaussian disk model to the PIONIER 
H-band observations, but an exact dynamical modeling using 
AMBER HR data and the HDUST code is postponed to the dedi- 
cated paper in preparation. A toy model that is unable to describe 
the dynamics of disk density waves would only be an unneces- 
sary or even confusing half-step and is not included in this paper. 



2. Observations 

2.1. Optical and near-infrared spectroscopy 

To identify long-term variations and trends in the spectroscopic 



intersects the disk. The large stellar v sin / of 
phes that the disk is equatorial. 



-400 km s im- properties, data obtained with the Herds spectrograph (Kaufer 



1998) between 1995 and 2003 that are included in the statistical 



Pott et al. (2010j) used the spectro-astrometric mode of the 



Keck interferometer to study the structure of the 48 Lib circum- 
stellar disk. The derived disk size was 1.65 mas in the K-band 
continuum, estimated at about seven times the stellar diameter 
The disk was also resolved in several Pfund lines. Although the 
spectral resolution did not exceed 1 000, the self -phase referenc- 
ing mode permitted the authors to obtain differential astrometric 
precision of 3 micro-arcsec. Pott et al. concluded that the contin- 
uum and the Pfund emission lines originate from a more com- 
pact region within the Bry emission zone, consistent with the 
fact that the Pfund lines are weaker than the Bry line. They also 
inferred that deviations from axisymmetry are radius-dependent, 
which is characteristic of one-armed disk oscillations. 

One purpose of this report is to put the recently resumed 
emission line variability in the circumstellar disk of 48 Lib 
into a common perspective with the general spectroscopic and 
polarimetric variability in the interval 1995-2011. We also 
present results of the first AMBER high-spectral resolution (HR, 
R - 12 000) interferometry of this star with a fairly complex 



study by Rivinius et al. ( 2006) were re-analyzed. This echelle 
spectrograph covers the spectral region of 3500-8600 A in 2 
arms with a gap of about 200 A between them. During the obser- 
vations of 48 Lib, the spectrograph was attached to the ESO (La 
Silla) 50-cm and Ondfejov 2-m telescopes. 

The ESO Science Archive only contains two sets of UVES 
( jPekker et al. 2000| l echelle spectra, which were obtained on 
July 11, 2001 and April 3, 2008. Although each of them is lim- 
ited to one night, they enable a detailed comparison of spec- 
tral lines before and after the density wave became detectable 
again (see Fig.[2|. Additional red and blue optical spectra at = 
16 000 were secured at the Pico dos Dias Observatory (OPD) of 
the Laboratorio Nacional de Astroffsica (LNA), Brazil, on July 
2, 2009 and May 12, 2011. The observations were made with 
a Cassegrain spectrograph equipped with a 1200 grooves/mm 
grating blazed at 6562 A and a 1024x1024 pixel CCD. To iden- 
tify long-term Ha V/R trends in 2007-201 1, we also used spectra 
from the BESS database (http://basebe.obspm.fr), which mainly 
comprises observations by amateur astronomers. The spectral 
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Fig.l. u-v coverage of the 48 Lib PIONIER and AMBER ob- 
sei-vations. Spatial frequencies corresponding to the minimum, 
central, and maximum wavelength are plotted for each baseline. 



resolution of some BESS datasets is quite low but they define 
the V/R variations well. 

In August 2009, a high-resolution (R - 45 000) spectrum 
of the Bry profile was obtai ned with the GEMINI/PHOENIX 
echelle near-IR spectrograph ( Hinkle et aLj200 3 ). The spectrum 



covers the interval 2.16-2.17 jum, and the 40 s integration gave a 
signal-to-noise ratio S /N ~ 170. 

The characteristics of the spectra and their origin are sum- 
marized in Table [T] 



2.2. Imaging polarimetry 

New multi-band linear polarimetry was obtained in the period 
May 2009-September 2010 using the lAGPOL imaging po- 
larimeter attached to the Cassegrain focus of the 0.6-m Boiler 
& Chivens telescope at OPD/LNA. We used the CCD camera 
with a polarimetric module described in Magalhaes et al. ( 1996 ), 
employing a rotating half-wave plate and a calcite prism placed 
in the telescope beam. A typical observation consists of eight 
consecutive wave-plate positions separated by 22°5. Details of 
the data reduction can be found in Magalhaes et al. ( 1984[ l. In 
each observing run at least one of HD 183143, HD 145502, and 
HD 187929 was observed as the polarized standard star to cali- 
brate the position angle. 

For comparison with the polarimetric status before the new 
V/R cycle commenced we used seasonal means in the U,B,V,R,I 
bands in 1990-1998 pubUshed by |McDavid| ( p00T||T994l[T999 i. 
The values are based on observations obtained with the po- 
larimeter described by Breger ( 1979), which is mounted on the 
University of Texas 0.9-m telescope at McDonald Observatory. 
For the period 1989-1994, they were complemented by V-band 
Halfwave Polarimeter data (HPOL; [Nordsieck & Harris|ri996| l 
from the 0.9-m telescope of Pine Bluff Observatory, which were 
retrieved from the web site http://www.sal.wisc.edu/HPOL 

All polarimetric observations were corrected for the inter- 
stellar polarization of (0.819±0.050)% and PA = (96 + 5)° de- 
rived by Draper & Wisniewski (private communication) using 
HPOL spectropolarimetry and a Serkowsky law with /Imax - 
5591 A. 



to 65 m. The fiinges were slightly dispersed over three spectral 
channels across the H band. The observing conditions were poor 
(DIMM (S arazin & Roddier 1990) seeing -1.5-2.0" and coher- 
ence time <i ms). But we were able to reduce the data with the 
standard pndrs package described by LeBouquin et al. ( 201 1\ . 
The observations provide the simultaneous measurement of six 
absolutely calibrated visibilities and four closure-phases in the 
H-band. 

High-spectral resolution (HR; R = 12000) interferometric 
observations in the K band were obtained with the AMBER in- 
strument (|Petrov et al. 2007) and the 8-m UT telescopes of the 
VLTI on May 12, 2009. For comparison, AMBER observations 
of ^ Tau obtained during science verification of the HR mode on 
October 7, 2008 were retrieved from the ESO Science Archive. 
The three-telescope baseline configuration was UT1-UT2-UT4 
for both stars. The observing conditions were average to good for 
48 Lib (DIMM seeing ~ 0.8-0.9", coherence time -6.3 ms), but 
quite poor during the ^Tau observations (DIMM seeing ~1.2- 
1.5", coherence time only ~ 1.5 ms). 

The central wavelength was set to 2. 171958 /im in the first 
order to match the Bry line. The Finito fringe tracker ( |Le| 
Bou quin et al.|2007 1 enabled us to record the whole spectral re- 
gion with integration times of 3 s and 1 s for 48 Lib and ^ Tau, re- 
spectively. HD 143033 (diameter 1.36 ± 0.04 mas) was observed 
as the calibrator after 48 Lib and HD 33833 (0.818±0.056 mas) 
before ^Tau. The diameters of both calibrators were adopted 



from the catalog by Merand et al. (2005 1. The dates of the ob- 
servations are given in Tab. |2| u-v plots of the 48 Lib PIONIER 
and AMBER observations are shown in Fig.[T] 

The reduction of the AMBER data was performed using the 
yorick amdlib3.0 package (TatuUi et al.|2007[|Chelli et al.|2009| l. 
The wavelength scale - affected by a problem with the position- 
ing of the AMBER grating - was corrected by fitting telluric 
lines using the Pyth on procedure writ ten by A. Merand. The 



same methods as in Stefl et al. (20091 were applied to derive 
absolute and diiJerential quantities for visibilities, phases, and 
closure phases across the Bry line. Following the recommenda- 
tion by Tatull i et al.] (|2007) and' Millour et al.| ('2008), only the 
top 50% of the frames with the best S /N were averaged for cal- 
ibrators as well as science targets to derive the transfer function 
and obtain absolutely calibrated spectra. However, because of 
the setting inconsistency between the science target and the cal- 
ibrator, only the differential visibilities and phases were derived 
for ( Tau. 

Differential visibilities, phases, and closure phase were aver- 
aged, and the values were normalized relative to the continuum 
in the intervals 2.150-2.158 //m and 2. 178-2. 190 jum. The differ- 
ential method enables a more precise measurement of quantities 
with respect to the adjacent continuum. Compared to the abso- 
lute method, a significantly higher S /N can be obtained because 
the atmosphere affects all spectral channels in nearly the same 
way. However, the absolute calibration of the angular size of the 
observed object is lost. 



2.3. Interferometry 

Broad-band interferometric observations in the H band were ob- 
tained with the PIONIER ( LeBouquin et al.|201 l| l visitor instru- 
ment and the four 1.8-m auxiliary telescopes of ESO's VLTI 
system ( GUndemann et al.|2003| l on June 04, 201 1 . The AT con- 
figuration was DO-Gl-HO-Il, with baselines ranging from 40 m 



3. Spectroscopic results 

3.1. Long-term V/R variations of emission lines 

The V/R ratio in our spectra was measured as the ratio of the 
emission peak heights (above the zero flux level) derived by fit- 
ting a Gaussian to the upper part of the peaks. 
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Fig. 2. a) Histoiy of 48 Lib Ha VI R variations (defined as ratio 
of the peak intensities above the continuum) in the period 1950- 
201L The circles mark original data from this paper and from 
Rivinius et al.|(2006[), + data compiled fro mfUnderhill (1953 ), 
Faraggiana (1969l) , |Gray & Mar lborough' (1974), Slettebak &, 
Reynolds ( 19 78 ), [Andrillat & Fehrenbach (1982) Hanuschik] 
etal. ( 1995 ) and |Banerjee et al | ( |2000| l."b)Hg and Bry VI R vari- 



ations (defined as the ratio of the full peak intensities) during 
the period 1995-2010. The arrows mark the dates of AMBER 
and UVES observations used for this study. Data from dif- 
ferent datasets are identified by the following symbols. Ha: 
+ [Rivinius et aL] ( [20061 HEROS spectrograph), x VLT/UVES, 
o BeSS database, □ LNA/Ecass; Bry: A VLTl/AMBER, ▲ 
Gemini/PHOENIX. c) Radial velocities of selected Hei lines; 
mean of Hei4026 and 5878 (o) and He 1 6678 (+). The RV er- 
ror is 3-5 km s"' for shell lines, but as large as 10-20 km s"' for 
He I lines due to their peculiar line profiles, d) Radial velocities 
of Fell (o multiplet 96, □ multiplet 117) and Sin, multiplet 2 
(+). 



Some VIR studies (e.g. [Okazaki[[T997[ ) report that 48 Lib 
went through four VIR cycles with lengths of about 13.5, 12, 
7.5, and 10 years between 1950 and 1994. The statement is sup- 
ported mainly by Fig. 1 in Hanuschik et al. ( 1995[ l, in which 
one can identify VIR maxima close to 1950 and approximately 
in 1963, 1975, 1982 and 1992. However, the plot combines VIR 
variations derived from Ha, Wfi and Hy lines and for the 1975 
maximum no Ha values are plotted. As indicated also by the 



large diff'erences between Ha and HyS/Hy values, combining VIR 
values of different Balmer lines may result in a misleading pic- 
ture both for observational and theoretical reasons: 

- Hy6 and y are formed in the inner disk closer to the central 
star than Ha and map different parts of the density wave spi- 
ral structure. Consequently, significant phase shifts between 
VIR variations of different Hi lines can b e expected - see 



Stefl et al.| ( [20091 ) for ^Tau or jKogure & Suzuki[ ( |1984| l for 



EWLac. 

- Intensities of the V and R peaks in H/8 and y are much more 
affected by a strong stellar absorption. No correction using 
synthetic or template profiles was performed in the early ob- 
servations. 

For the above reasons a recompilation of the VIR variations in 
Ha VIR only was performed and combined with the new 1995- 
2011 data. Results are shown in Fig. l2^. Because for some of 
the published datasets the authors used VIR defined as VIR* - 
{I(V) - I(cont))l(I{R) - Ucont)) and did not pubHsh Ha profiles, 
this definition had to be accepted also in Fig. l2^. Unlike in Fig. 
Igl^, \ogiyiR*) is plotted. 

The plot reliably, although not precisely in time, shows max- 
ima in 2005 and 1992. A weak maximum in 1980/1981 agrees 
with that in Hanuschik's plot. No Ha observations were found 
in the literature for the period 1963-1976. Fig. l2^ shows a few 
higher values around Hanuschik's maxima in 1950 and 1963, but 
also a large scatter All VIR values before 1985 were measured 
from published Ha profiles derived from photographic plates. 
We conclude that only two Ha VIR cycles are sufficiently docu- 
mented in the star's history, giving cycle lengths of about 1 1 and 
15 years. Only the last cycle, see Fig. I2l5,c,d, will be analyzed 
in more detail in this paper 

Unfortunately, even for the last cycle, there are gaps just 
when the quasi-periodic oscillations were probably absent 
(1995-1998) and re-started (at some time between 2004 and 
2007). Although detailed information about the Ha behavior in 
these phases is, therefore, missing, 48 Lib offers one of the not 
so frequent opportunities to study the recovery of VIR activity in 
a Be star 

Lower amplitudes seem to be connected with a shorter cy- 
cle length. The observations obtained in 2000 show a first small 
increase in VIR and thereafter almost constant VIR ~ 1.0 for 
more than two years. A strong increase to VIR » 1 must have 
taken place at some time between 2004 and 2007. Beginning 
with the first observation on JD54 172 after this gap, VIR con- 
tinued to decrease and now reached the 1996 level of about 0.8. 
This places a lower limit on the length of the descending branch 
of the peak in the VIR curve of about 4.2 years. Taking into ac- 
count the history of the variations since 1970, see Fig. l2^, we 
can assume that the variability in Fig. ^ is cyclic. From the 
speculative assumption of a symmetric curve, a VIR cycle length 
of about 17 years could be derived. This length would be sig- 
nificantly longer than the previous ones observed in the second 
half of the 20th century. The present descending VIR branch also 
suggests - consistently with the longer cycle - that the descent 
started from at least the same or even a higher value than the 
maximum in the previous activity cycle. 

The upper left panel in Fig. |3l shows that the H/? VIR ratio 
follows a similar trend as the one of Ha. 

Apart from Ha VIR variations. Fig. l^l^ also includes Bry 
VIR measurements made in spectra extracted from the AMBER 
data and from the single PHOENIX HR observation. If the Bry 
VIR curve exactly matches that of Ha, one could infer that Bry 
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is ahead of Ha by about 300-330 days. For a mean cycle length 
of 15 years, this would correspond to a phase shift of 19-22°. 



3.2. Absorption line profiles 

Shell absorption lines in Be stars are formed in the narrow sec- 
tor of the circumstellar disk projected on the central star. Since 
the involved range in azimuth is much larger along the stel- 
lar equator than in the more extended circumstellar disk, which 
moreover - in the Keplerian case - rotates much more slowly 
than the star, circumstellar absorption lines are mostly much 
narrower than photospheric lines. In the simplest case of an 
axi-symmetric, Keplerian disk, the line profiles should be time- 
invariant. 

The shell-star nature of 48 Lib manifests itself through nu- 
merous strong and relatively sharp lines from singly-ionized 
metals in its spectrum. Given the above basic model of Be-shell 
stars, their temporal variability is evidence of large deviations in 
the disk from a circular Keplerian radial velocity (RV) field with 
Vrad = Okms"^ or a strong azimuthal asymmetry in the density 
distribution, or both. Their profiles are an important indicator of 
the disk dynamics, which they sample in a complementary way 
to emission lines. Furthermore, Hanuschi k & Vrancken| ( |1996| l 
discovered extremely narrow multiple components of Fe n shell 
lines variable both in RV and depth on a timescale of less than 
one day. The authors hypothesized that these features could be 
caused by co-orbiting clumps in the disk or by higher-order dis- 
tortions of the global density wave causing the V/R variations 
seen in the emission lines. 

Fig. [3] shows line profiles of Hfi and some shell and He i 
lines in selected epochs in 1995-2011. The shell lines exhibit 
variations in depth and width, which are linked to V/R (previous 
subsection) and RV variations (see next subsection) in Ha. Both 
depths and widths of the shell lines reach a maximum close to 
the phase of minimal V/R and RV. Owing to missing data, the 
observational correlation is not clear for V/R > 1. The deeper 
lines may also show a blue-wing asymmetry as it would result 
from motions toward the observer. It is worth noting that Fig. [3] 
may be somewhat misleading because it is biased by the lower 
resolution of the OPD/LNA spectra. 

The profiles of He i lines are quite unusual. The entire optical 
region covered by the UVES and HEROS spectra contains only 
two lines with profiles not obviously strongly affected by su- 
perimposed shell components, namely He 1 47 13 and He 1 5786. 
However, even the core of the latter is probably contaminated 
by circumstellar absorption. In spite of the low S/N in He 1 47 13, 
Fig. [3] illustrates that the line suffers from strong profile and 
asymmetry variations. 

He 1 4471 and most other He i lines show fairly complex pro- 
files with one or more absorption components. Most of them turn 
out to be superimposed metallic lines, e.g. Fen4469.485 and 
4472.086, or Fen 6677.305 in the lines He i 4471 and 6678, re- 
spectively. However, because most He i lines follow the same RV 
trend as the shell lines, one cannot exclude that their absorption 
cores are blends of metal and He i lines. Extended blue wings 
of He 1 5876 and 6678 can be recognized during MJD52 000- 
54 000, while depressions in the red wing occurred in the interval 
MJD 50 000-5 1 000 and around MJD 55 000. 

The observed line-profile variability strongly suggests that 
the above assumed distinction between photospheric and cir- 
cumstellar lines is not nearly as simple as presented. On the other 
hand, the profiles of lines formed in the disk are important line- 



of-sight probes of the disk dynamics, which they sample in a 
complementary way to emission lines. The same diagnostics is 
not available for Be stars not viewed through the disk. 

The only datasets of this study with spectral resolution suf- 
ficient for the detection of very narrow (FWHM~ lOkms ') 



optical absorption components as reported by Hanuschik & 
|Vrancken| ( | 199 6) are those from UVES. However, no such nar- 
row component could be detected in Fe ii profiles. A search in the 
HEROS spectra, the very first of which was obtained only two 
days after those by Hanuschik & Vrancken on March 10, 1995 
(their Fig. 1) also remained without positive results. Because the 
spectral resolution of HEROS at 5300 A is only ~15 km s \ one 
cannot decide whether the components were just not detected or 
really disappeared. A splitting of the absorption core or traces of 
extra absorption components can be seen in some Fe ii lines in 
the HEROS spectra obtained in 1996, but they cannot be reUably 
attributed to the same phenomenon. 

3.3. Radial velocities 

Variable line profiles are quantitatively easiest characterized by 
radial-velocity (RV) measurements, especially if they have nar- 
row cores. Without higher-order perturbations, RVs of photo- 
spheric and shell lines in Be stars should be the same. However, 
|Rivinius et al | ( |2006| l pointed out very high negative radial ve- 
locities of shell lines (-200km s ') in 48 Lib although the star 
probably is not a binary (cf. Sect|6]l. Timescales of years are dif- 
ficult to reconcile with any stellar processes. Therefore, the vari- 
able radial velocities result from gas motions in the disk, and, 
as shown by Okazaki (1997), the cycle lengths match periods of 
disk oscillations well. But there is no straightforward interpre- 
tation of the RV of line cores in terms of particle trajectories or 
streaming motions because the line profiles are convolutions of 
velocity profiles with density and excitation profiles. 

To study the nature of the time-dependent RV field, radial 
velocities of Feii and Siii lines were measured by fitting a 
Gaussian to the cores of their absorption profiles. The results 
plotted in Fig. |2|l show that Hff V/R and shell RV variations are 
well correlated. The stagnation of the V/R variations in 1996- 
2003 was accompanied by an increase in the RV of shell lines 
by about 100 km s"'. In 2011, the values returned to the same 
negative velocity of about -100 kms"' as in the end of the 1990s. 
Separately plotted in Fig. |2jl are mean velocities for Feii mul- 
tiplets 96 (spectral lines at 4508, 4549 and 4583 A) and 1 17 
(lines 5197, 5234, 5275, and 5317 A) and Sin multiplet 2 (fines 
6347, 6371 A). In the LNA spectra, only lines of Feii multiplet 
96 and Si ii could be measured. In spite of this inhomogeneity 
and different spectral resolutions, the values for different multi- 
plets and for Feii and Sin do not show significant differences. 
Accordingly, the differentiation of the velocity field with excita- 
tion potentials may be small. 

The contamination of He i profiles by metallic lines implies 
the potential risk that the He i RVs are biased. The risk was min- 
imized by 

- Checking for the presence of metallic lines in the region of 
each He I line using a catalog of spectral lines (Reader et al. 
|1980| l. He 1 6678, formed in the upper photosphere or even 
inner disk, as well as He 1 4026 and He 1 5876 were identified 
as the least contaminated ones according to both line catalog 
and visual inspection of the lines. The He i4713 line was not 
included in the mean. This line is very shallow so that RV 
measurements are strongly affected by noise and imperfectly 
corrected echelle ripple functions. 
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Fig. 3. Long-term variations of Hy8 and selected Feii, Si ii and He: line profiles obtained in 1995-2011. Note the very different scale 
of the normalized flux, indicated at the bottom of each panel. All HEROS and UVES spectra were smoothed with a 5-pixel filter 
The spectrograph and date are marked in the right margin. 



Measuring He i RVs in the line wings by fitting the profile be- 
tween approximately one and two thirds of the core intensity 
by a Gaussian. The FWHM of the shell lines (variable, 50- 
80kms ' as measured in HR UVES spectra) is much lower 
than that of Hei lines (300-380 km s '). Provided that shefl 
lines do not aff'ect the line wings, they bring about a line 
profile distortion and a larger scatter in RV, but they do not 
significantly affect the measured RVs. The fact that the RVs 
of different He i lines are consistent also excludes the pos- 
sibility that they are mistaken for those of shell lines. The 
corresponding shifts of the entire broad He i absorption pro- 
files can also be recognized visually in Fig. |3] 



The errors are largest for the two lines that are supposedly 
mainly formed in the photosphere. He 1 4026 shows two com- 
ponents, which are much broader than dynamical shell lines. 
The blue component dominated the RV measurements in 1995- 
2003 and was included in the means plotted in Fig. |2] Similar, 
but much more difficult to measure, components are present in 
He 1 4009. 

Apart from the large scatter, there are also considerable sys- 
tematic diff'erences between the RVs of individual lines. Note 
that the values plotted in Fig. [2] are not representative of the 
radial velocity of He i. They were merely included to illustrate 
some basic similarities to conventional sheU lines. 



Somewhat surprisingly, the RVs of selected Hei lines broadly 
agree with those of more classical shell lines. Fig. [2]: shows 
RVs of He 1 6678 and the mean RV of Hei 4026 and Hei 5876. 
The scatter in the RVs of the He I lines is considerably higher 
than for shell lines. Fig. [3] also documents that line cores of all 
He I lines are variable on a long timescale exhibiting wide ranges 
in shape and strength. All these effects lead to errors in the RVs 
of He I lines of a few tens of km s ' . 



4. Imaging polarimetry 

[McDavidl ( |1999[ l performed a time-series analysis of his broad- 
band polarimetry covering more than a decade. While the polar- 
ization angle, 0, was constant, the degree of the polarization, p, 
suggested quasi-cyclic variability, which McDavid fitted with a 
period equal to one-half of the cycle length of the Ho- VjR cycle. 
The peak-to-peak amplitude in p was about 0.1%. A correlation 
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with the VI R variability could be caused by the density varia- 
tions associated with one-armed disk oscillations. The double- 
frequency variability of the polarization might arise because po- 
larization cannot distinguish between the density enhancement 
being on the left or the right side of the disk. Unfortunately, the 
OPD/LNA polarimetry does not cover the VjR stillstand, which 
might have provided an interesting test of the proposed, still un- 
confirmed interpretation. 

The upper panel of Fig. |4] shows a decreasing V and / 
band polarization trend in 2009-2010. A similar decrease is 
present also in the remaining spectral bands. The relative er- 
ror of an individual observation is 2-A%, while the drop of po- 
larization in the given period is 10-15%. The lower two pan- 
els compare the mean values of p and for 48 Lib in the pe- 
riod 1990-1998 derived from the data of ' McDavid| ( fT994l \\^9[ 
with the OPD/LNA observations obtained in 2009-2010. V- 
band Pine Bluff Observatory data from 1989-1994 yield slightly 
lower values of p than those derived from the McDavid dataset: 
p(V) = (0.519 + 0.070)% and 0(y) = (150.3 ± 1.0)°. Because for 
an optically thin or moderately thick disk, the polarization angle 
is always perpendicular to the disk plane ([ Wood et al.||1996| , 
the value of derived by McDavid is equivalent to the disk po- 
sition angle (60.3±1.0)°. The errors show standard deviations, 
which are biased by real intrinsic variations as shown in the up- 
per panel. Similarly, the 1990-1998 means are derived from sea- 
sonal averages and may be affected by variations on a timescale 
of years. 



Assuming overestimated standard deviations due to real 
trends, the difference between the mean polarization in the 90s 
and 2009-2010 can be considered as real and together with the 
decreasing temporal trend suggests that the polarization was in- 
creased at the time of the V/R recovery and in 2009-2010 it was 
returning to the level measured in the 90s. To within the errors, 
there is no change in the wavelength dependence. The maximum 
polarization appears in the B band. 



5. Interferometry 

5.1. Continuum 

The disk continuum emission around 48 Lib is marginally re- 
solved by the H broad-band PIONIER interferometric observa- 
tions. We studied its morphology via parameter modeling us- 
ing the LITpro package ( Tallon-B osc et al. 2008 ). Results are 
summarized in Tab. [3] To facilitate the comparison with other 
published disk diameters, we first fitted the dataset with a uni- 
form disk model, which yields ^ud = ( 1 .34 ± 0.05) mas. We then 
modeled more realistically the H band emission as a Gaussian 
profile disk and an unresolved component (the central star) that 
contributes two thirds of the total flux ( [Dougherty et al.|[l994| l. 
The derived FWH M of 6gd - (L4 4 + 0.05) mas is consistent 
with the results by [Pott et al. ( 2010| l (1.65 mas) obtained in the 
K-band continuum with the Keck interferometer 

Although an axisymmetric Gaussian disk plus an unresolved 
component reproduce the full PIONIER dataset relatively well 
(X^ = 2.07), we explored the possibility to constrain the disk 
position angle. We obtained a better fit (x^ = 1.33) with an 
elongated Gaussian disk whose major axis has a FWHM of 
= (1.72 + 0.2) mas and a position angle of (50 ± 9)° (north 
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Fig. 4. 48 Lib polarization observations: Upper panel - decreas- 
ing p trend in V (o) and I (+) bands in the period 2009-2010. 
Lower panels - mean polarization and polarization angle vs. 
spectral band in the periods 1990-1998 (black A) and 2009- 
2010 (red o). Straight lines in the bottom panel show the regres- 
sion fits to the BVRI bands for both epochs. Within the errors 
the spectral slope does not change. 
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polarimetric and the interferometric PA of the disk are consis- 
tent. 

Finally, we explored the possibility of a fully resolved com- 
ponent that would be the signature of a dim shell with a spatial 
extent significantly larger than the disk itself. The resulting best- 
fit model has a negligible shell contribution of 0.01 + 0.02 of the 
total flux and an increased value compared to the model with- 
out shell. We conclude that a shell whose integrated contribution 
is more than a few percent in the H band is not probable around 
48 Lib. 
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Fig. 5. Absolute visibilities for 48 Lib in the H band obtained 
with PIONIER and auxiUary telescopes at the DO-HO-Gl-Il 
VLTI baseline configuration. Black dots mark the observed vis- 
ibilities, crosses (x, red) the fit with a Gaussian disk, and plus 
signs (+, gi^een) the fit with an elongated disk. 
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Fig. 6. Relative photocenter shifts derived from the AMBER 
differential phases across Bry. After excluding the outliers at 
extreme offsets and points corresponding to the continuum (- 
50<NE<50) the linear regression fit (blue line) to the photocen- 
ters across the line (plotted in blue) yields the disk position angle 
(48.5+1.5)°. The NE shifted points correspond to the shorter- 
wavelength wing of the Bry emission profile, the SW shifted 
points to its longer-wavelength wing. 



5.2. Bry line 

The spectro-interferometric observations of 48 Lib were carried 
out with AMBER during the descending branch of the V/R ac- 
tivity (see arrows in Fig. |2]). Because Fig. |2]does not cover the 
48 Lib V/R maximum and the cycle length is much more vari- 
able than in f Tau ( jStefl et al .12009) 1, the phase estimate is little 
accurate. One can, however, estimate that the AMBER observa- 



Table 3. Best fit for the parametric modeling of the H-band con- 
tinuum emission from PIONIER data. Position angles are de- 
fined north to east. 



Model 


Parameter 


Value 




Xr 


Unifomi Disk 




1.34 H 


I 0.05 mas 


2.17 


Star + Gaussian Disk 




(1.44 


± 0.05) mas 


2.07 


Star + Elongated 




(1.72 


± 0.20) mas 


1.33 


Gaussian Disk 


•^GD / ^GD 


(1.66 


±0.3) mas 






p^major 


(50 ± 


9)= 




Star + Elongated 


^GB 


(1.55 


± 0.20) mas 


1.35 


Gaussian Disk 


^major //iminor 
GD ' GD 


1.93 H 


= 0.60 




+ dim shell 


p^major 


(53 + 


10)° 






excess/totalflux 


0.01 H 


= 0.02 





tions took place close to phase 90° since V^R. This is also the 
approximate phase of the AMBER observations of ^ Tau, which 
justifies the inclusion of the latter data for comparison. 



The high-quality observations of 48 Lib permit good-fidelity 
differential visibility and phase profiles to be derived (Fig. |7]i. 
The disk was well-resolved at all three baselines, both in contin- 
uum and line emission. At the UT1-UT4 baseline, the disk was 
almost over-resolved in Bry. Although the data were absolutely 
calibrated, the analysis of continuum visibilities yields results 
inconsistent with PIONIER continuum visibilities and with po- 
larimetric observations. Consequently, only differential visibili- 
ties and phase profiles over Bry are discussed below. 

For disk emission lines, simple uniform or Gaussian disk 
models predict a single-peak central reversal (akin to a spectral- 
line flux profile with superimposed circumstellar emission line) 
in the visibility and a smooth S-shaped phase profile. Indeed, 
Fig. |7] baseline UT1-UT2, exhibits such an S-shape phase pro- 
file. However, overall the 48 Lib HR AMBER profiles are sig- 
nificantly different. All Bry visibility profiles show a double- 
component structure with a central peak and an extra peak, su- 
perimposed to the S-shape phase profile conventionally seen at 
lower spectral resolution. 

The height of the central peak in the 48 Lib visibility function 
reaches the level of the continuum at the intermediate baseline 
and about half of it at the longest and shortest one. The central 
visibility peak seen at the UT1-UT4 baseline is accompanied by 
a multi-peak phase profile. At the UT2-UT4 baseline, only a very 
weak phase reversal emerged in the negative lobe of the S, and 
the UT1-UT2 baseline shows a pure S-shaped profile (see Fig. 

The low-noise 48 Lib data also reveal a rich structure in the 
closure phase profile. The high amplitude of the closure phase 
indicates a strong asymmetry. However, the splitting into both 
positive and negative peaks is well above the noise level. 

In spite of the higher noise, a similar double-component 
structure can be seen in the ^Tau visibility and phase profiles 
(Fig. [sj. The central peak in almost reaches the continuum 
level at the UT2-UT4 baseline. Unlike in 48 Lib, the visibility 
shows a gradual increase in the longer-wavelength wing and an 
asymmetry between the two visibility components. For ^Tau, 
both interferometry and polarimetry give a disk major axis po- 
sition angle of 32°( |Stefl et al.||2009| ). This is almost parallel to 
the UT1-UT4 baseline and the difference in position angle to the 
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Fig. 7. Differential visibilities, differential phases, and closure phase for 48 Lib. Lengths and position angles of individual baseUnes 
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Fig. 8. Differential visibilities, differential phases, and closure phase for ^Tau. The baselines were the same as for 48 Lib (cf. Fig. 
0. 



other two baselines is only about 30°. Although none of the base- 
lines was close to perpendicular in the October 2008 AMBER 
observations. Fig. [Hjstill indicates some features over and above 
the basic S-shape. But their significance should be verified by 
better quality data with different baseline angles. 

Not only ai-e the HR AMBER observations of 48 Lib and 
^Tau quite similar, they are also fully consistent with compa- 
rable observatio ns of five other classical Be stars discussed by 



Stefl et al.| ( |20T2] i. 



The differential phase signal can be converted into the 2D as- 
trometric shift piA), which defines shifts of the photocenters in 
individual wavelength bins across the spectral line and the ori- 
entation of the disk major axis projected on the sky. Here we 
consider the applicability of the photocenter approach to 48 Lib 
and other Be stars. The photocenter can be interpreted using the 



Zernike-van Cittern theorem and moments of the flux distribu- 
tion ( Lachaumel|2003 see their Eq. 4 for the definition of the 
nth moment). In this formalism, the phase of the interfero met- 
ric signal is proportional to the first flux moment, defined as 
the photocenter vector, and a third-order term, related to the ob- 
ject symmetry. The symmetry is described by the closure phase 
(Lachaume's Eqs. 8b and 10). If the closure phase is close to 
zero, the conversion from the phase to the astrometric shift can 
be made by inverting the formula 



cf>(A)^-2n-p(A)-^ 

A 



(1) 



where b is the interferometric baseline vector projected onto the 
sky. Because the phase does not depend on the second moment in 
the first approximation a sufficient condition for the conversion 
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is a negligible closure phase, not necessarily a marginal visibil- 
ity, as required for the mean diameter or symmetry expressions. 

A density wave in the Be star disk obviously leads to asym- 
metries and moderate or even large closure phase, see Fig. |7]for 
48 Lib. Consequently, Eq. 1 can be used only as a rough estimate 
of the disk orientation. 

Fig. |6] shows the relative photocenter offsets across the Bry 
line with respect to the continuum photocenter. Although it over- 
constrains the problem, the astrometric shifts were computed for 
all three baselines and their consistency was checked. The 2-D 
photocenter was reconstructed from the two shortest baselines 
(UTl-2, UT2-4) with the lowest noise and sufficiently different 
PA, see Fig. [T]. The UTl-4 baseline is inconsistent. The rela- 
tively good agreement of the disk position angles derived from 
the PIONIER fit and the AMBER photocenter shift indicates that 
the shorter baselines in AMBER observations were still close to 
the astrometric regime. The approaching part of the rotating disk 
(the blue wing in the Bry line) is oriented in the NE direction, 
the receding part to the SW. 

6. Discussion 

Observations of 48 Lib between 1995 and 20 11 document a re- 
covery of the VIR activity and suggest a coiTelation between the 
underlying density wave and RVs and profiles of metallic shell 
and He I lines. During a long plateau in the VIR ratio through 
~2003 the RV still continued to increase so that an absence of 
detectable VjR variations does not exclude ongoing global disk 
oscillations. 

The one-armed disk oscillation model has so far not been 
used to simulate a case like the one observed in 48 Lib with 
substantial asymmetries of shell lines apparently caused by very 
large shifts of their cores. Such high RV amplitudes of shell lines 
were also observed in other Be-shell stars, e.g., ^Tau ( Riviniusl 
let al.|2006| l. Nevertheless, in the model of Carciofi et al. (2009F 
for the Be star ^ Tau, the particles orbit the star in very eccen- 
tric elliptical orbits. In that particular model, the radial veloci- 
ties are as high as 100 km s"'. Elliptical orbits could also explain 
the strongly variable width of metallic shell lines. Obviously, 
such large intrinsic variations are much more readily observed 
in disks with an edge-on orientation towards the observer, i.e. in 
Be-shell stars. 

48 Lib is also of special interest because many of its He i 
lines also exhibit RVs in line wings, which broadly track those of 
shell lines. Because circumstellar components of He i lines must 
form very close to the central star, they are a very useful diag- 
nostic of the inner structure of the disk. Observations of Be stars 
( Rivinius et al. 1998[ l and modeling of their light curves ( |Carciofi 



|et al.||2011| l provide increasing evidence for a time-dependent 
mass decretion rate and an episodic refilling of the disk. After 
the decretion is stopped or significantly reduced, the time scale 
to stabilize the inner disk in the resulting low-density state is 
only on the order of weeks to months. The observed strengths 
and velocities of the He i line cores suggest that a high density is 
maintained in the inner disk through a high decretion rate. The 
large outer disk radius estimated from interferometry would be 
in line with this interpretation. 

Extended blue wings of He 1 5876 and 6678 as observed 
during MJD 52 000-54 000 could also be a mass-loss indicator 
However, the observed dependency on VjR phase, i.e., azimuth 
angle of the disk, makes this less likely. An alternative inter- 
pretation is offered from our preliminary modeling of the disk 
dynamics. In some VIR phases, the line-of-sight radial velocity 
shows an extra minimum in the very inner disk. Such a velocity 



field could then be responsible for the extended blue wings and 
even doubl e shell profil es, which were observed close ioV - R 
in f Tau by Stefl et al. ( |2009) [their Fig. 3 and related discus- 
sion]. 

In Be double stars, orbital and VIR periods are in some cases 
known to be synchronized (Stefl et al._2007j . Their orbital pe- 
riods are on the order of tens to hundreds of days, e.g. 4 Her: 
46.18 d ( |Koubsky et al.|1997l l or eCap: 128.5 d ( iRrvmius et al. 
2006). For a putative binarity of 48 Lib, Fig. |2] would be con- 
sistent with a period of 15-25 years. The shape of the RV curve 
would imply only minor eccentricity. Therefore, a binary would 
stay far away from the disk at all times and the very extreme disk 
variability could hardly result from the presence of a companion 
star, making a binary model very unattractive for 48 Lib. 

Three independent methods - polarimetric angle, LlTpro 
modeling of continuum visibilities and astrometric analysis of 
the photocenter shifts - yield the position angle of the major disk 
axis in the interval 48°to 60°. It is worth mentioning that each of 
the methods may suffer from some systematic errors. The polari- 
metric angle of 48 Lib depends strongly on the correction for the 
interstellar polarization. The elongated axially symmetric disk 
assumed in the LlTpro solution can be considered to be only the 
first-order approximation. The astrometric analysis was applied 
under conditions that do not guarantee the convertability from 
the interferometric phase of the used AMBER observations and 
must be considered only as a rough estimate. The three methods 
provide information on parts of the disk at different distances 
from the central star While the polarization comes from the in- 
ner disk region, which is also close to the region of formation of 
the H-band continuum, the Bry line is formed much farther out 
in the central part of the disk. Considering all the above facts, the 
data do not exclude the possibility that the position angles of the 
inner disk and the H-band continuum and Bry-emitting regions 
are nearly identical. 

The accuracy of the AMBER absolute visibilities reaches 
about 5% under optimal conditions (Chelli et al. 2009) and can 
be estimated between 6 and 10% in our 48 Lib observations. The 
inconsistency of the AMBER absolutely calibrated continuum 
visibilities with those from PIONIER, polarimetric observations 
as well as simple LlTpro models indicate other problems with 
the absolute calibration of the AMBER HR continuum visibil- 
ities. Possible, but so far little investigated sources of system- 
atic errors are the fringe tracker FINITO and integration times as 
long as several seconds. Note that the test by Chelli et al. ( 2009| l 
was made only for AMBER medium resolution and integration 
times between 0.02 and 0.16 s. 

A reliable calibration of the continuum visibilities could be 
made for the PIONIER H continuum visibilities. Their fitting 
by the relatively simple LlTpro disk models confirms the disk 
diameter of 1 .7 mas determined by the Keck H broad-band in- 
terferometry and the position angle of 60°derived from our po- 
larimetric observations. This diameter at the Hipparcos parallax 
of 6.97 mas corresponds to 0.25 AU. Assuming a B3V main- 
sequence star radius of 3.56 7?© ( Ha rmanec|| 1 998| l, the H con- 
tinuum disk diameter is about 15/?*. According to Fig. 1 in 
ICarciofi (2011 1, the K-band continuum is formed typically at 
about 1-3/?*, but Bry only about three times farther from the 
central star. Consequently, we can assume that the diameter of 
the Bry emitting region might be at least 30-50/?*. 

The LlTpro estimate confirms that the 48 Lib cicumstellar 
disk has one of the largest angular diameters among Be stars. 
A consistent modeling of the 48 Lib disk using more compre- 
hensive interferometry and spectral energy distributions is un- 
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der consideration. It would also test the Keplerian dynamics as- 
sumed by the disk oscillation model, for which the extreme ra- 
dial velocity amplitudes of shell lines may still be a challenge. 

Only thanks to the AMBER high-spectral resolution mode 
could the double-component structure of visibility and phase 
profiles be discovered. In 48 Lib, the Bry line width is about 
45 A so that atR = 12 000 about 25 wavelength bins sample the 
central visibility maximum very well and can reveal even more 
subtle features in the visibility profile. At AMBER's medium 
spectral resolution (R ~ 1500) just three wavelength bins can- 
not establish the same structure. In spite of the very high phase 
stability, no double peak structure or traces of a central peak can 
be seen in visibility profiles of |Pott et al. (2010 their Fig. 2), 
obtained about a year before the AMBER data and atR = 2000. 
Spectral resolution is also of key importance for interferometric 
studies of the disk dynamics. 

The occurrence of similar visibility and phase pr ofiles in the 
more y/7?-stable disk of ^Tau as well as other stars (Stefl et al. 
2012) 1 shows that they do not result from a dynamical peculiarity 
related to the V/R recovery in 48 Lib. Instead, these structures are 
a general property independent of the phase of disk evolution. 

Comparison with the 48 Lib Bry emission profile offers a 
simple qualitative explanation of the W-shape visibility profiles. 
The Bry visibility is much reduced in the line center, where the 
line emission drops and the visibility is given mainly by the con- 



tinuum visibility. jKraus et al. ( 201 1^ analyzed HR interferomet- 
ric observations of the Be stars ^Tau and ySCMi and proposed 
that phase reversals correspond to the transition from the first to 
the second visibility lobe. This over-resolution effect can explain 
the phase reversals in the case of large and close disks like that of 
48 Lib. This interpretation is consistent also with the fact that the 
phase reversal at the UT2-UT4 baseline and the second closure 
phase minimum appear at the position of the first and deeper vis- 
ibility minimum. However, the lobe transition may not provide a 
general explanation of th e phase reversal, p articularly for less re- 
solved and distant disks. Stefl et al. (2012jl speculated that phase 
reversals and multi-peak closure phase profiles can be caused 
by secondary effects in the disk and may hint at additional mor- 
phological disk features. However, until a realistic modeling of 
the disks including the density waves is performed, the general 
explanation of the HR phase and closure phase profiles remains 
open. 



7. Conclusions 

The compiled spectroscopic data of 48 Lib confirm that its cir- 
cumstellar disk is evolving on a timescale of years to decades. 
The correlation between the Ha V/R ratio and RVs of the cores 
of metallic shell and even He i lines implies a dense inner disk 
fairly close to the photosphere or even touching it. An elongated 
disk model fit of the broad-band interferometric data indicates an 
outer disk radius as large as 15 ^^in the H continuum. This leads 
to an estimate of the Br y radius of several tens of stellar radii. 
Together with a relatively large parallax of 6.97 mas, these re- 
sults make 48 Lib a prime candidate for spectro-interferometric 
studies of the dynamical structure of Be-star disks. A more exact 
model would also shed light on the question whether the high- 
velocity cores of shell and He i lines find their explanation within 
the density wave model alone or if an episodic gas outflow to, or 
infall from, the circumstellar disk must be included. 

The resumption of the V/R oscillations in the first decade of 
this century was accompanied by a deepening of the disk shell 
profiles and an increased polarization degree. 



For more than 20 years, the polarization angle has been con- 
stant to within the errors and remained unaffected by the disk 
oscillation, which, therefore, must be purely planar without ma- 
jor 3-D precession or warping. 

At high spectral resolution, interferometric visibility and 
phase show significant deviations from the simple S-shape pro- 
file conventionally adopted for the interpretation of interfero- 
metric data of Be stars at low and medium spectral resolution. 
Only a realistic modeling (Faes et al., in preparation) will be 
able to answer the question if double-component phase or clo- 
sure phase profiles are always only a consequence of interfero- 
metric effects during visibility lobe transitions or if they (also) 
reflect secondary effects in the disk. 
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